Arp 302: Non-starburst Luminous Infrared Galaxies 

K.Y. Lo, Yu Gao, and Robert A. Gruendl 
Laboratory for Astronomical Imaging, 
Department of Astronomy, University of Illinois, 
1002 West Green Street, Urbana, IL 61801 



Received 



accepted 



- 2 - 



ABSTRACT 

Arp 302, a luminous infrared source (Lir = 4.2 x 10 11 L ), consisting of two 
spiral galaxies (VV340A and VV340B) with nuclear separation of ~ 40", has 
the highest CO luminosity known. Observations with the BIMA array at 5" 
x7" resolution reveal that the CO emission is extended over 23.0 kpc in the 
edge-on spiral galaxy, VV340A, corresponding to 6.7 x 10 10 M of H 2 . In the 
companion face-on galaxy, VV340B, the CO emission is extended over ~ 10.0 
kpc, with 1.1 x 10 10 M Q of H 2 . The large CO extent is in strong contrast to 
starburst systems, such as Arp 220, in which the CO extent is typically < 1 
kpc. Furthermore, L IR /M(H 2 ) is found to be £ 6.0 L /M throughout both 
galaxies. Thus the high IR luminosity of Arp 302 is apparently not due to 
starbursts in the nuclear regions, but is due to its unusually large amount of 
molecular gas forming stars at a rate similar to giant molecular clouds in the 
Milky Way disk. Arp 302 consists of a pair of very gas-rich spiral galaxies that 
may be interacting and in a phase before a likely onset of starbursts. 

Subject headings: galaxies: individual (Arp 302, VV340, UGC 09618) - galaxies: 
interactions - galaxies: spiral - galaxies: starburst - infrared: galaxies - ISM: 
molecules 
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1. INTRODUCTION 

Ultra-luminous infrared (IR) galaxies (L IR > 10 12 L Q ) and most luminous IR galaxies 
(LIRGs; 10 12 L > L IR > 10 11 L ) are almost always found in interacting or merging 
systems (e.g., Sanders et al. 1988; Melnick & Mirabel 1990; Murphy et al. 1996). On the 
other hand, only a small fraction of optically selected interacting/merging galaxy systems 
are truly IR-luminous (e.g., Arp 1966; Arp & Madore 1977; Kennicutt et al. 1987; Bushouse 
et al. 1988). Apparently, galaxy interaction is a common causal factor for the enhanced IR 
luminosity, but it is not the only one (cf. Lawrence et al. 1989; Leech et al. 1994). 

As the large IR luminosity can be explained by a burst of star formation and almost 
all LIRGs are found to be molecular gas-rich (Sanders et al. 1991; Solomon et al. 1996), 
understanding of LIRG mergers then depends on how interaction induces starbursts and 
enhances the star formation efficiency in these merging systems. A study of the gas 
properties at various phases of the merging process would help identify the key physical 
processes involved. For example, all four ultra-luminous IR galaxies and most LIRGs in 
advanced merging stages observed with sufficiently high spatial resolution, show a high 
CO surface brightness in their nuclear region with an extent <^ 1 kpc typically. This has 
been interpreted as a very high central H 2 surface density, usually more than 10 4 M pc~ 2 
(Scoville et al. 1991; Yun & Scoville 1995). 

Few studies have been able to explore the gas-rich galaxy mergers over various merging 
phases (e.g., Combes et al. 1994). In a study of 50 LIRG mergers, the CO luminosity 
is found to increase with increasing separation of the merging nuclei, which may imply 
that the molecular gas content of LIRGs is decreasing as merging advances (Gao 1996; 
Gao & Solomon 1996). Previous studies have concentrated on relatively more advanced 
merger systems in which the interstellar medium (ISM) has already been highly disrupted 
by the starbursts. In order to isolate the conditions in the ISM leading to starbursts, we 
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are starting a program to study a sample of LIRGs chosen to represent different phases 
of the interacting/merging process, using the newly expanded Berkeley-Illinois-Maryland 
Association (BIMA) millimeter-wave array (Welch et al. f 996) which is ideally suited to 
this problem given its large primary beam and wide spectral bandwidth. The goal is to 
sample statistically the evolution of the properties of the molecular clouds in the galaxies 
as they go through the merging processes. 

Here, we present the first results of this program on Arp 302. Arp 302 (UGC 09618) 
consists of two galaxies (VV340A and VV340B - Sb and Sc, with apparent magnitudes 
14. 6 m and 15. 3 m , respectively) with a separation of ~ 40" between the two nuclei. At a 
median redshift (cz) of 10160 kms -1 , its luminosity distance, dL, is 136.6 Mpc (H = 75 
kms _1 Mpc _1 , q = 1/2) and it has an infrared luminosity, L IR , of 4.2 x 10 11 L . 

Single-dish measurements show a CO luminosity L C o ~ 1-9 x 10 10 K kms -1 pc 2 , which 
corresponds to a molecular gas mass of ~ 9 x 10 10 M , the largest amount of molecular gas 
known among the samples of ~ 50 LIRG mergers studied recently by Gao (1996) and ~ 37 
ultra-luminous galaxies surveyed by Solomon et al. (1996). Given the large IR luminosity 
and the close proximity of the two galaxies, the normal assumption would be that Arp 302 
is powered by starbursts that result from the interaction of the two galaxies. The relative 
large separation (~ 25 kpc = ~ 40" at the angular distance of 136(l+z)~ 2 = 128 Mpc) 
compared to merger remnants such as Arp 220 would further suggest that Arp 302 is in an 
early phase of merging/interacting. 

2. OBSERVATIONS AND DATA REDUCTION 

Observations of Arp 302 were started with the BIMA 6-element array in October 1995. 
More observations were made with the new 9-element array in the C configuration in April 



- 5 - 



and B configuration in May 1996 (Welch et al. 1996). The digital correlator was used 
with the widest bandwidth to cover a velocity range of 1,600 kms -1 , at a resolution of 8.4 
kms" 1 . The ~ 2' primary beam was centered between the two galaxies to cover all possible 
emission in this interacting system which extends slightly more than 1'. The observing 
frequency was 111.4952 GHz, corresponding to the CO (1-0) line at the median redshift 
cz=10,160 kms -1 of this merging galaxy system. 

The nearby quasar 3C345, used for the complex antenna gain calibration, was observed 
for 8 minutes before every 30 minutes on Arp 302. The strong quasar 3C273 was observed 
at the beginning for 30 minutes to calibrate the spectral dependence of the gain. The 
typical system temperature ranged between 300 and 800 K. 

The interferometer data were reduced with the MIRIAD data reduction package (Sault 
et al. 1995). Observations in May 1996 on baselines longer than 50 kA were found to 
suffer from atmospheric decorrelation and were therefore weighted by applying a Gaussian 
taper (effective resolution of 3."5) in the UV plane prior to inverting the UV data cube. 
Data cubes of the CO brightness distribution, T b (x, y, V), where x and y are the spatial 
coordinates and V is the radial velocity, at various velocity resolutions (8.4, 20.0 and 50.0 
kms -1 ) were created covering the range from —600 to 600 kms -1 about cz = 10,160 kms -1 . 
The cleaned maps have a final synthesized beam of 7" x 5", with the major axis of the 
beam at a position angle of 10°. The various moment maps were then made to produce the 
integrated intensity map (Ico( x )y) — JT b (x, y, V)dV in K-kms" 1 ), the intensity weighted 
velocity map, and the intensity weighted velocity dispersion map. 



3. RESULTS 
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3.1. Spatial Distribution of the CO Emission 

In Figure 1 (Plate), the contours of an integrated CO intensity map at 7" x 5" resolution 
are superposed on a R-band CCD image of Arp 302 (left panel); whereas the right panels 
show the CO spectra from various locations in Arp 302. The Ico distribution is very well 
correlated with the R-band CCD image of the galaxies, especially along the prominent dust 
lane in the northern edge-on galaxy VV340A. The CO emission is especially bright in the 
inner ~ 14", or ~8.7 kpc. At the angular distance of 128 Mpc, 1" corresponds to 620 pc. 

About ~ 85% of the total CO flux (S co = (2k/ A 2 ) / T b (x, y, V)dxdydV in Jy kms" 1 ) 
is from the edge-on galaxy, distributed along the major- axis over ~ 35" and marginally 
resolved along the minor axis. If the CO emission is coplanar with the galactic disk, the 
shape of the CO distribution would place a limit on the inclination angle, i, of the galaxy 
to be >75°. 

While ~ 15% of the total CO flux is from the southern galaxy VV340B, the CO surface 
brightness is weaker, well resolved by the 7" x 5" beam, and appears correlated with the 
central bar. Some correspondence of weak CO emission with the northern spiral arm is 
also indicated (Fig. 1). Within the uncertainties, all the CO flux detected from these two 
galaxies with the NRAO 12m telescope (Gao 1996) is detected in the BIMA maps. 

We have also detected CO emission from a position away from the two galaxies at the 
position: a(J2000.0) = 14 /l 57 m 01. s 6, 5(J2000.0)=24°36'48" (cf. spectrum labelled Off in Fig. 
1) that does not correspond to any optical feature. The reality of this interesting source has 
to be verified by further observations. 
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3.2. Kinematics 

Figure 2 shows a position-velocity diagram of the CO surface brightness in VV340A 
along the major axis (PA ~ 5°) with ~6" and 20 kms -1 resolution. It shows that the 
kinematics of the molecular gas in VV340A has largely a constant velocity gradient (85.5 
kms -1 kpc" 1 ) within the central 14", or ~8.7 kpc, along the major axis and a velocity 
range of 780 kms -1 at FWZI. The faint CO emission to the north of the central 15" (~ 10 
kpc) appears to have a constant velocity between 260 kms -1 and 344 kms" 1 , and to the 
south between —180 kms" 1 and —210 kms" 1 . At the extreme southern end of the disk, 
the integrated intensity map (Fig. 1) shows a slight deviation westward from the linear 
morphology exhibited by the rest of the galaxy. In the same region, the velocity field also 
has a distortion from the rotational signature. This could be interpreted as evidence that 
the interaction between the two galaxies has perturbed the molecular gas distribution. 

A constant velocity gradient has often been described as indicating a "solid-body" 
rotating disk. In the case of the inner 14" of VV340A, the "solid-body" disk would have 
a 4.35 kpc radius. Alternately, a 8.7 kpc diameter rotating molecular ring, while also 
consistent with the constant velocity gradient, would be easier to understand dynamically. 
A ring of molecular gas may also explain the shape of the CO spectrum from VV340 A 
(Fig. 1), and the weakness of CO emission near the systemic velocity (~ 40 kms _1 in the 
figures, or cz = 10200 kms" 1 ) from the nuclear region. Thus, a strong concentration of 
molecular gas in a ring structure rotating with a circular speed of 390 kms" 1 (corresponding 
to the inner disk) combined with a thin, extended, molecular gas disk with a flat rotation 
velocity of ~240 kms" 1 viewed edge on would be consistent with the observations. We 
note, however, that the observed systemic velocities of the rotating ring and disk, if taken 
as the centroid of the extent of each respective component, differ by ~30 kms -1 . This may 
arise if the distribution of material in the plane of VV340A is not axi-symmetric and has 
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non-circular motions. Further study of the gas kinematics is necessary. 

There is a velocity gradient across VV340B as well, although the total velocity extent 
across VV340B is only ~ 100 kms _1 (FWZI). Given the large amount of molecular gas 
and the narrow line-width, the galaxy must be nearly face-on, which is consistent with the 
shape of the galaxy in the R-band image (Fig. 1). 

From the centroid of the CO line profiles for VV340A and B, we determine the relative 
systemic velocity between the two galaxies, Va — Vb = —340 ± 20 kms -1 , with VV340B 
blueshifted relative to VV340A. Thus, the spin vector of VV340A could be nearly parallel 
to the orbital angular momentum vector of VV340B, unless a very large relative transverse 
motion exists. 



3.3. Molecular Mass and Surface Density 

Arp 302 is the most CO luminous system known in the local universe, more luminous 
in CO than those ultra-luminous galaxies out to the redshift of ~ 0.27 (Solomon 
et al. 1996). The total molecular gas mass of each galaxy can be estimated from the 
total CO flux, using M(H 2 ) = 1.18 x 10 4 dL 2 Sco, where dL is the luminosity distance 
in Mpc and Sco is the total CO flux in Jy kms -1 . We use the conversion factor, 
X = [N(H 2 )/Ico] — 3-0 x 10 20 cm _2 /(Kkms _1 ) for convenient comparison to previous 
results, even though it is not a priori applicable to Arp 302. 

The H 2 gas mass in VV340A is (6.7 ±0.3) x 10 10 M , with >85% in the central 8.7 kpc. 
The dynamical mass interior to the rotating molecular ring (4.35 kpc radius) is estimated 
to be ~ 1.5 x 10 11 M Q . This implies that if He is included, the gas mass constitutes 
^60% of the dynamical mass of VV340A within a radius of 4.35 kpc. The H 2 gas mass in 
VV340B is (1.1 ± 0.1) x 10 10 M Q . The combined H 2 gas mass in this merging system is 



- 9- 



(7.9 ±0.3) x 10 10 M . 

The HI mass in Arp 302 has been measured to be ^,2.3 x 10 10 M (Mirabel & Sanders 
1988). Thus, the total gas (HI + H 2 + He) content of Arp 302 is - 1.4 x 10 11 M . Even if 
X is 5 times smaller, the total gas content would still be ~ 0.54 x 10 11 M . However, this 
does point to the importance of a direct determination of the conversion factor X for Arp 
302. 

Although the southern end of the VV340A molecular disk is distorted, no molecular 
gas has been detected between the two galaxies in the disk overlap region. But, the CO 
source to the east of VV340A has a molecular mass of (1.0 ± 0.2) x 10 9 M . 

The spatial extent of the CO source in Arp 302 provides estimates of the face-on H 2 
surface density, E H2 . Thus, we obtain E Ha of ~ 10 3 M pc~ 2 within the radius of 4.35 kpc 
and ~ 20 M pc~ 2 in the outer disk in VV340A, and ~ 100 M pc~ 2 for the inner part of 
VV340B. In the inner parts of both galaxies, S H2 is high compared to that of the Milky 
Way disk (4 — 17 M pc~ 2 ; Sanders et al. 1984), but small compared to that of the LIRGs 
and ultra- luminous galaxies at the more advanced merger stage (typically > 10 4 M pc~ 2 ; 
e.g., Yun & Scoville 1995). 

4. DISCUSSION 

4.1. Is Arp 302 a Star-burst System ? 

While the L m = 4.2 x 10 11 L certainly qualifies Arp 302 as a LIRG, the CO emission 
is extended over the entire ~23 kpc diameter disk in VV340A and over many kpc in 
VV340B. The large CO extent is in stark contrast to most of the more than 20 LIRGs and 
ultra-luminous galaxies for which high spatial resolution CO (1-0) observations have been 
made. They showed highly concentrated CO emission in nuclear regions <^ 1 kpc in size 
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typically (Scoville et al. 1991). 

Spatially unresolved observations of LIRGs showed that they tend to have very large 
Lir and Leo, usually interpreted as indicating very high star formation rate and H 2 gas 
mass, as well as very high star formation efficiency (i.e., L IR /M(H 2 ); Solomon & Sage 1988; 
Young et al. 1989; Sanders et al. 1991; Solomon et al. 1996). From spatially resolved 
observations, one can derive local measures of CO surface brightness both in terms of T b 
and I C o converted to S Ha , the H 2 surface density. For LIRGs and ultra-luminous galaxies, 
all these quantities are extremely high compared to Galactic values (e.g., Scoville et al. 
1991; Sargent & Scoville 1991; Yun et al. 1994; Yun & Scoville 1995) 

The ratio of Lir/M(H 2 ), used to indicate the star formation efficiency, ranges between 
20 to 100 L /M for nuclear starburst regions (Scoville et al. 1991). For Arp 302, the 
global ratio is 5 L /M . If we scale the IR luminosity and extent with those of the radio 
continuum emission from Arp 302 (Condon et al. 1990; Condon et al. 1991), we obtain 
values of 6, 2.5 and 2 L /M for L IR /M(H 2 ) in the inner 8.7 kpc region and the outer part 
of VV340A, and VV340B, respectively. This is comparable to the ratio for giant molecular 
clouds (GMCs) in the Milky Way (L IR /M(H 2 ) ~ 4, Scoville & Good 1989) but much below 
that observed in Arp 220 (95 L /M ) and other starburst systems (cf. Scoville et al. 
1991). Furthermore, we can infer a mean IR surface brightness (Sir) of 7000, 55 and 300 
L pc -2 for the same three regions respectively. This is very low compared to the typical 
~ 10 5 L pc -2 in nuclear starburst regions (Lo et al. 1987, and cf. Scoville et al. 1991). 

Therefore, the high IR luminosity in Arp 302 arises simply because an extremely large 
amount of molecular gas is forming high mass stars at a rate similar to the GMCs in the 
disk of the Milky Way and not because starbursts are occurring in the galactic nuclei. 
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4.2. Is Arp 302 an Interacting System ? 

The two galaxies (VV340A and VV340B) in Arp 302 have a projected nuclear 
separation of 40" and a velocity difference of 340 kms -1 . Are the two galaxies physically 
associated, and not just chance superposition ? It has been shown via Monte Carlo 
simulations that the posterior probability of two galaxies brighter than 15 m with a 
velocity difference of between 200 and 1000 kms™ 1 being bound is 0.88 (Schweizer 1987). 
Furthermore, the distortion of the CO distribution in the southern end of VV340A may be 
an indication that the two galaxies are in fact interacting. The HI emission line profile from 
Arp 302 (U 09618 in Mirabel & Sanders 1988) is significantly different from the CO profile, 
in that HI emission appears to be missing at the velocities corresponding to VV340A. 
Future VLA observations of HI in Arp 302 may show direct evidence for the interaction of 
the two galaxies in terms of HI tidal tails and the kinematic signature. 

We can also estimate the ratio of the kinetic and potential energies of the interacting 
system: T/U = (M r V 2 /2)/(GM A M B /R), where M r is the reduced mass, V is the relative 
velocity, and R is the nuclear separation of the two galaxies with masses, Ma and M B . If 
M A w M B as implied from their comparable extent, M r ps M b /2, and T/U = RV 2 /(4GM A ). 
For Ma = RaV 2 /G, where V r is the rotational velocity of VV340A at the outermost radius 
R A , T/U = (V/V r ) 2 (R/R A )/4. Since V r « 250 kms' 1 with V >340 kms" 1 (if the transverse 
motion is taken into account) and R ^ 2 x R A (the projected separation), it is plausible 
that T/U is on the order of unity and E = T — U ~ 0. (If the mass of one galaxy is 
substantially less than the other, the T/U ratio would change by a factor of two.) This 
analysis suggests that the two galaxies could be on parabolic orbits about each other. 

The relative radial systemic velocity of the two galaxies is 340 kms -1 and the spin of 
VV340A is parallel to the apparent orbital motion. This indicates that the disk of VV340A 
may be undergoing a prograde encounter with its perturber VV340B, unless there is a very 
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large transverse motion between the two galaxies. Numerical simulations have shown that 
two gas rich galaxies of comparable mass released on parabolic orbits, with the mass and 
orientations similar to those found in Arp 302, would merge within two orbital times, or a 
few times 10 8 years (Barnes and Hernquist 1991; see also White 1979; Toomre & Toomre 
1972). 

Thus, we conclude that the high probability that the two galaxies in Arp 302 are 
bound, the distortion in the CO distribution in VV340A, the plausibly small total energy 
of the two galaxies, and the numerical simulation results are all indications to suggest that 
Arp 302 is an interacting system. 

Given both the stars and the gas distribution in Arp 302 system do not appear 
highly disturbed, the two galaxies are apparently approaching each other for the first time, 
suggesting that Arp 302 is in an early phase of the interaction. The next phase of Arp 302 
may be represented by the galaxy VV114 (Yun et al. 1994), in which the molecular gas has 
merged into a common envelop before turning into an ultra-luminous galaxy. 

5. SUMMARY 

With the BIMA array, we have resolved the CO emission from the luminous IR galaxy 
system Arp 302 into galactic scale distributions. Local measures in Arp 302, such as 
Sh 2 , L ir /M(H 2 ), Sir, are all much smaller than those found in starburst systems. Thus, 
the enormous IR luminosity of Arp 302 is apparently not due to starbursts, but is due to its 
unusually large amount of molecular gas forming stars at a rate similar to giant molecular 
clouds in the Milky Way disk. 

Also, evidence suggests that the two galaxies in Arp 302 are interacting and at an 
early phase of the interaction. As the interaction progresses, bursts of star formation may 
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eventually be initiated. If so, we may be witnessing two interacting gas-rich galaxies in a 
pre-starburst phase. Identifying and studying such systems are important for isolating the 
conditions in the ISM that will lead to starbursts. 
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Fig. 1. — Integrated CO(l-O) intensity (Ico) contours are overlaid on a R-band image 
(grayscale) of Arp 302 (VV340A and VV340B). The contours shown are 6, 8, 12, 24, 48, 
and 96 x 3.2 K-kms" 1 with a 7" x 5" beam. To the right of the image, the top and bottom 
spectra, derived from the data cube, show the flux density (S„ in Jy) of the CO emission 
from the two entire galaxies VV304A and VV304B, respectively. The four spectra shown 
between the top and bottom spectra are the flux density of the CO emission from a single 
beam area centered as follows: North (a=14 /l 57 m 0. s 8, 5=24°37 / 20"); South (a=14 h 57 m 0. i 55, 
5=24°36'48"); Off (a=14 ft 57 m l. s 65, 5=24°36'48") ; and Arm (a=14 /l 57 m 0. s 2, <5=24°36'34"). 
kms -1 corresponds to cz = 10,160 kms -1 . 



Fig. 2. — Posit ion- Velocity diagrams for the two galaxies VV304A (a) and VV304B (b) 
in Arp 302 system. The contours for both panels represent T b , starting at 0.064 K and 
increasing by steps of 0.064 K. The first white contour is at 0.51 K. 
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